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Abstract 

We suggest a theory of internal coherent tunnehng in the pseudogap region, when the apphed voltage U is 
below the free electron gap 2Ao. We address quasi ID systems, where the gap is originated by spontaneous 
lattice distortions of the Incommensurate Charge Density Wave (ICDW) type. Results can be adjusted also 
to quasi-lD superconductors. The instanton approach allows to calculate the interchain tunneling current 
both in single electron (amplitude solitons, i.e. spinous) and bi-electron (phase slips) channels. Transition 
rates are governed by a dissipative dynamics originated by emission of gapless phase excitations in the course 
of the instanton process. We find that the single-electron tunneling is allowed down to the true pair-breaking 
threshold at Uc = 2Was < 2Ao, where Was — S/ttAq is the amplitude soliton energy. Most importantly, the 
bi-electronic tunneling stretches down to Uc — (in the ID regime). In both cases, the threshold behavior is 
given by power laws J ^ {U — Uc)^ , where the exponent (3 ^ vp /u is large as the ratio of the Fermi velocity 
vp and the phase one u. In the 2D or 3D ordered phases, at temperature T < Tc, the one-electron tunneling 
current does not vanish at the threshold anymore, but saturates above it at U — Uc ^ Tc Aq. Also the 
bi-particle channel acquires a finite threshold Uc — W2tt ^ Tc Aq at the energy of the 2tt phase soliton. 
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I. INTRODUCTION. 



A. Pseudogaps and the subgap tunneling. 

Interchain, interplane transport of electrons in low dimensional (quasi ID, 2D) materials attracts 
much attention in view of striking differences between longitudinal and transverse transport mech- 
anisms, revealing a general problematics of strongly correlated electronic systems Beyond a 
low field (linear) conduction, the whole tunneling current-voltage J-U characteristic J{U) and the 
conductivity a = dJ/dU are of particular importance. A common interest in tunneling phenomena 
in quasi ID conductors has been further endorsed by recent experimental achievements, see and 
especially in new techniques of the intrinsic tunneling (see [6] for a short review) where electronic 
interchain transitions take place in the bulk of the unperturbed material. Particularly appealingjs 
the access to topologically nontrivial excitations - solitons, showing up in subgap spectra [3, [s 



Role of preexisting solitons in tunneling spectra of High-Tc materials has been advocated in 
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loi. 



The goal of this article is to describe t 



Charge Density Wave (ICDW) , see ^, Isl, 1^ . The following content of this chapter will introduce 
experimental evidences, theoretical grounds and conditions for existence of solitons. Chapter II 
will describe the theoretical approach. Chapter III will give results of calculations, which details 
can be found in the Appendix. 

If the tunneling processes were going between free electron states, as they have been formed by 
the rigid ICDW, then the current onset would correspond to the voltage = 2Ao of the gap in the 
spectrum of electrons. But actually there is also a possibility for tunneling within the subgap region 
Eg < U < Eg. It is related to the pseudogap (PC) phenomenon known for strongly correlated 
electrons in general, well pronounced in quasi ID systems, arid particularly in cases where the gap 
is opened by a spontaneous symmetry breaking (see 



leir manifestations and peculiarities in Incommensurate 
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|l5[ and refs. therein). The 



PC is originated by a difference, sometimes qualitative, between short living excitations which are 
close to free electrons, and dressed stationary excitations of the whole correlated systems. The 
true excitations gain the energy, but loose the probability of transitions to these complex states. 
In the ICDW, the dressing of a bare single particle results in self-trapped states like the amplitude 
solitons with energies Was < below the free electron activation energy Aq; then the subgap 
interval 2Was < U < 2Aq will be observed as a PC Below this pair-breaking threshold 2Was 
(observed in optics, spin susceptibility), in tunneling there may be also contributions of collective 
states - phase solitons, which energy Wps is even much lower. These states are activated in the 
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bi-electronic channel - coherent tunnehng of two electrons. In the ID regime (which we shall mostly 
address below) Wps = 0, then the subgap tunneling fils the whole pseudogap. 

B. Solitons in ICDWs 

Strongly correlated electronic systems show various types of symmetry breaking originating de- 
generate ground states (GS). The degeneracy gives rise to topologically nontrivial perturbations 
exploring the possibility of traveling among different allowed GSs. Of special interest are totally lo- 
calized and truly microscopic objects, namely solitons (or instantons for related transient processes) 
which can carry single electronic quantum numbers: either spin 1/2 or charge e, or both. Being 
energetically favorable with respect to electrons, solitons would determine the electronic properties 
which can be proved theoretically at least for quasi ID systems. We address here the cases of a 
continuous GS degeneracy of an Incommensurate CDW (ICDW) (see ) where solitons can 

be created in single items via single electronic processes. Results can be generalized also to quasi- 
ID superconductors. Continuous degeneracy of the ICDW order parameter Acos{Qx + if) comes 
from an arbitrary chosen phase 99, which characterizes the freedom of translation 6x = —5ip/Q of 
the ICDW as a whole. At a first sight, the CDW is commonly viewed as a narrow gap anisotropic 
semiconductor with most of its properties described by free electrons e or holes h near the gap edges 
ibAo. Thus, Aq would give the activation energy in kinetics and thermodynamics characteristics 
(conductivity, spin susceptibility, heat capacitance, NMR); the same Aq would be observed from 
dynamic probes as photoemission, external tunneling. The double gap 2Ao would be measured as 
the edge in optics or in internal tunneling described below. Doping or junction injection (FET) 
would require the threshold Aq and lead to formation of electronic pockets near the gap edges. But 
actually almost nothing from this well established semiconducting picture takes place in ICDWs. 

1. Activation energies (Ay and A^) measured from transport in the on-chain (Ay) and in the 
interchain (A 1 ) directions differ by several times (typical values for TaSa are Ay ~ 200 K and 
A^ ~ 800 K [16]) which signifies an intrinsically different character of processes involved; 

2. Energies deduced from the spin activation A^ = 640 K 13] and from the dynamic relaxation 
Ad Q are in between ^ 600 K; 

3. Optical absorption shows the peak at the scale of 2A_l - thus associated to 2Ao, but its shape 
is strongly smeared out with intensity deeply spread within the expected spectral gap [l^ - the 
pseudogap effect; 

4. Thresholds for charge transfer may not exist at all, or have low values compatible with Ay - i.e. 
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associated with the interchain decoupHng scale, see 
5. Charge injection results in absorption of additional electrons into the extended GS via so-called 
phase slip processes (see rather than in formation of Fermi pockets. The static phase slip - a 



2tt soliton as an on-chain CDW defect, has been directly observed in recent STM experiments 




FIG. 1: Soliton trajectories in the complex plane of the order parameter. The circle is the manifold of 
allowed ground states |A| — Aq. The vertical diameter line 20 = vr is the stable amplitude soliton. The 
vertical chordus line is an intermediate soliton within a chiral angle 26. The value 26 — 100° is chosen, 
which corresponds to the optimal configuration for the interchain tunneling (see [isi l and below). The arcs 
with arrows show phase tails required to level out the perturbation at large distances from the soliton. 

The order parameter of the ICDW is the complex field A = A{x,t) e'x.p[iip{x,t)], acting upon 
electrons by mixing states near the Fermi momenta points zizkp. Solitons in the ICDW rise from 
strong interactions between electronic e, h and collective (amplitude A and phase ip) degrees of 
reedom, which is set up since all of them are modified by spontaneous formation of the new GS 



22l |. Electrons' dressing by the collective deformations results in relatively heavy self-trapped 
particles with lower energies A < Aq, which originates two basic scales. The high one can be 
observed, as 2Ao or Aq, for instantaneous processes like optics, tunneling, photoemission when the 
ground state has no time to adapt itself to a perturbation. The lower activation A can be observed 
for fully relaxed dressed particles, which might dominate stationary effects of thermodynamics. 



kinetics, NMR, etc., see refs. in [12] • The GS degeneracy 
elementary excitations: they become topological solitons 



gives a special character to these complex 



22( 1 ■ These are the trajectories connecting 
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different equivalent GSs; tliey are cliaracterized by a cliiral angle 29 = ip{x = +00) — ip{x = —00). 
Tliis angle defines a family of "chordus solitons" - FiglU see Appendix for details; it provides a 
continuous path in the configurational space for electron's self-trapping. The time evolution 9{t) 
describes the self-trapping dynamics - the instanton [3], see FigJH] below. 




FIG. 2: Selftrapping branches (in units of Aq) of chordus solitons: total energies V„{6) for midgap fillings 
V = 1^2 oi localized split-off state and its energy level i?o = Aq cos(6') as functions of 9. 

Processes of self-trapping of electrons or their pairs are important in several respects. They 
provide: 

1. Spectral flow by transferring the split-off state Eq between the gap edges Aq — Aq, FiglJl 

2. Particle flow (conversion of normal carriers to the collective GS) between the conduction and 
the valence bands, which crosses the gap riding upon the split-off state; 

3. Microscopic Phase Slips, by adding/subtracting the 27r winding of the order parameter. 



FIG. 3: Profiles for the amplitude soliton: oscillating electronic density, its overlap, and the spin density 
distribution due to the unpaired electron at the midgap state. 



For single-particle processes (see reviews 



23j ) the self-trapping of one electron at Aq or 



of one hole at — Aq will proceed gaining the energy until the configuration takes a stable form 
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of the amplitude soliton AS: A{x) = Aq tanh(x/^o) cos((5x + ip), FigJSl Now £'0 = becomes 
a pure midgap state occupied by a single electron; thus the soliton carries the electronic spin 
s = 1/2. Curiously, the electric charge is zero rather than e, being compensated by the dilatation 
of electronic wave functions of the filled band E < — Aq. It can be interpreted in a way that 
the AS is symmetric with respect to the charge conjugation: it is the adaptation of an electron 
added to the GS of 2M particles, as well as of a hole upon the GS of 2M+2 particles. (Breaking of 
charge conjugation can be still treated within exact solutions [2^ - it gives a small electric charge 
to the AS.) Thus, in a ID system, the AS is a realization of a spinon, the particle carrying the 
elementary spin 1/2 but no charge. This particle plays a central role in phenomenological pictures 
of strongly interacting electronic systems in general l|. Theoretical value for the AS energy is 
Was = 2/7rAo ~ 0.65Ao; thus, the energy « O.SSAq is gained by converting the electron into the 
soliton. 

Excitation of the solitonic state takes place in those moments, when collective quantum fluctu- 
ations create an appropriate configuration with the necessary split-off intragap state [12]- At first 
sight, one needs to prepare the AS in its full form (the diameter in the Fig{T]), the probability of 
which is very low. But actually the spontaneous deformation is more shallow. Indeed, the energy 
Eg is yielded by the transfer of a single electron; hence the split-off energy level, prepared for 
this electron by the optimal fluctuations, must be at Eq = Was] therefore, the tunneling takes 
place when quantum fluctuations accumulate the chiral angle such that cos 9 = 2/tt, which gives 
29 « 100°, rather than 180°. 

The tunneling threshold 2Was of the pair-breaking is required only for the single particle channel, 
while for the bi-particle channel there may be no need for an activation as if the system is still 
metallic. This expectation resides upon the fact that the GS of the ICDW can incorporate a finite 
even amount of particles without paying an activation energy, unlike for odd numbers. For a finite 
length, the energy cost will be the same as for the metal: for each accommodated pair it is ^ hvp/L 
- vanishing for the large system length L. The accommodation requires for the total CDW phase 
increment along the chain Aip = 2tt per the added pair (see the branch V2 at the Figl2]). For the 
zero energy charge exchange between chains, it is necessary to adjust the number of states to the 
number of particles: that is at the chain a to transfer the emptied level at — Aq just below the 
gap across the gap up to the manyfold of properly empty levels above -I-Aq; the same time at the 
chain b the empty level at -|-Ao must be drawn down across the gap to form a new level at — Aq 
which will accommodate the arrived pair. This job is done by a simultaneous phase slip processes 
of opposite directions: negative ^ — 27r at the chain o, and positive ^ 27r at the chain b. With 
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the vanishing, in ID regime, of the phase shp energy Wps — > 0, the optimal fluctuation will require 
for the exactly mid-gap split-off state Eq 0, i.e. for 9 vr/2, which is the AS. Hence the AS 
becomes an optimal configuration, the barrier, for nucleation of the phase slip. 

Effects of 2D,3S long range order. In systems of correlated chains the local configuration 
must conform to the long range order. At long distances the phase ip must return to the mean 
value (ifoo = in our choice of Figl3|). 

The pure AS (vertical diameter) is accompanied by phase tails (arc lines) adjusting the phase 
discontinuity 7r/2 — > —tt/2 to the asymptotic zero value. With the total increment 2 x 7r/2 = vr of 
the phase concentrated within these tails, they carry the electric charge e, in addition to the spin 
1 /2 concentrated within the AS core. The composed particle carries both of the electron quantum 
numbers, but they are localized at very different scales. Contrary to the spin, the charge is not 
associated to a particular electron level; the electric charge transfer by the tails of the AS is rather 
the collective conduction promoted by ASs. The price for the adaptation to the long range order 
is an increase of the AS energy by an amount of the order of the transition temperature Tc- For 
the phase soliton it becomes the scale of its total energy Wps ~ Tc. 



C. Electrons and solitons in coherent tunneling experiments 



Recent tunneling experiments 



performed on quasi-lD materials with CDWs, namely NbSea 



and TaSs, can be interpreted as the first direct observation of solitons in dynamics. Tunneling 



spectra of NbSes 



show sharp peaks which are identified with the intergap tunneling at 2A2 and 
2A1 for the two CDWs coexisting in this material (see figures lb and 2 in [6]). Additional subgap 
features appear for both CDWs, when the zero bias contribution of normal electrons is suppressed 
by high magnetic field or by temperature T (see figure 2 in [Q]). These quite strong peaks scale 
together with A(T) at Vas ~ 2/3A. This unexpected feature finds a natural interpretation in the 
above picture of solitons as special elementary excitations. V^^ appears as the threshold for the 
tunneling between the Fermi level (middle of the gap) states of the electronic pocket, specific to 
NbSes, and the nonlinear multi-electronic complex - the Amplitude Soliton. 

The tunneling conductivity is also observed at even much lower voltages: U > Vt (see figures 



la,3a,5a in 



[sD, which can be provided only by the bi-electronic channel [l^ . 



In the ICDW it 



is realized via 27r phase solitons, which correspond to one period of stretching/squeezing of one 
chain with respect to the surrounding ones. These are the elementa ry p articles with the charge 
ib2e and the energy Wps ~ Tc (recall their recent visualization in (2lf]), which gives the right 
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order of magnitude both for the deep subgap tunnehng range and for the longitudinal conductivity 
activation. The tunneling process corresponds to a coincidence of oppositely directed phase slips at 
neighboring chains, which transfer a pair of particle between their GS. A theory of these processes 
will be presented below in the section IIIICI 



II. SELF-TRAPPING PROCESSES. 



A. Techniques. 



We shall follow the adiabatic method 
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13l . Q, ll^, developed for the X-ray absorption 



(PES, ARPES), optics, and for tunneling in low symmetry systems supporting polarons. It as- 
sumes a smallness of collective frequencies, here of 2Kf phonons ujph, in comparison with the 
electronic gap: ujph <^ Aq. Now, electrons are moving in a slowly varying potential A(x,t), so that 
at any instance t their energies Ej{t) and wave functions ^j{x,t) are defined from a stationary 
Schroedinger equation H^(x,t) = E{t)^{x, E(t)). The Hamiltonian H = H{x,A{x,t)) depends 
on the instantaneous configuration A{x,t), so that E{t) and ^(x, E{t)) depend on time only para- 
metrically. The time will be chosen along the imaginary axis where the necessary saddle points of 
the action are commonly believed to be found (see e.g. {29!). (Examples of profound analysis of 
the complex time contours for typical self-trapping problems can be found in the review 30|.) 

The process is determined by the change z^, = 0,1,2, of the number of electrons; it is a 
relative filling factor of the split-off intragap state (z^ > for an electron and < for a hole). At 
a given v, the instantaneous configuration A{x,t) determines the action S{i') = S{i',A). In the 
imaginary time, this functional is given as a sum of the kinetic and the potential Vu energies: 




S{u, A) = I dt{ I dx^^-^ + K(A(x, t)) \ -,¥^ = ¥0 + \u\E (1) 

dep^ph J 

where gep is the electron-phonon coupling constant and Uph is the bare phonon frequency. The 
potential contains the energy of deformations g~p J Ap, and the sum over electron energies in 
filled states (counted with respect to the GS energy for the total number 2M particles). It includes 
both the filled vacuum states at E < —Aq and the split-off intragap ones £"0, — Aq < Eq < Aq. 
In the non perturbed, A = Aq, state = |z^|Ao: the particle added to the non deformed GS, is 
placed at the lowest allowed energy, Aq. More details are given in the Appendix. 

The self-trapping evolution of the added electron can be fortunately described by the known 
(see [2_3;]) exact solution for intermediate configurations characterized by the singe intragap state 
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Eq = Aq cos 9 with < 8 < TT. It is known to be the Chordus Soliton with 29 being the total chiral 



angle: A(+oo)/A(— oo) = exp(2i6'), see Fig|2] and detai 
In calculations we follow the approximation 
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s in the Appendix. 



15( 1 of the zero dimensional reduction 



Aj{x,t) Asoiix - Xj{t),9j{t)) ; S{Aj{x,t)) (2) 

which limits the whole manyfold of functions Aj(x,t) to a particular class of a given function Agd 
of X (relative to a time dependent center of mass coordinate Xj and the local phase (pj). The 
function Ag^i is chosen as the chordus soliton trajectory which gives just one and only one split-off 
state Eq = Aq cos 6. 

The total on-chain action becomes 

S = S{ti,t2,U) = Score + Sdis (3) 

were Sdis is the dissipative contribution ([8]) described below and in the Appendix, while Score is 
the local action from the core of the chordus soliton 

Score = r dt[f{0)e^ + v,{e)\ (4) 

J tl 

Here 2f{9) is the variable "effective mass" of the "0-particle": 



f{e)= [ dx^V 



dA, 2 



de 



(5) 



and Vy{9) is the potential branch (see Appendix and Figl2|). Its basic form is 

v^{e) = AQ[{\u\--e)cose + -s\iie] (6) 

vr TT 

As the topologically nontrivial object, the chordus soliton cannot be created in a pure form 
with a finite chiral angle. In general, adaptational deformations must appear to compensate for 
the topological charge. These deformations are developing over long space-time scales and they 
can be described in terms of the gapless mode, the phase alone. Hence, allowing for the time 
evolution of the chiral angle 9 9{t) within the core, we should also unhinder the field (p <f{x, t) 
at all X and t. Starting from x — > — oo and finishing at x — > oo, the system follows closely the 
circle |A| = Aq, changing almost entirely by phase - notice the arc lines at the Fig(TJ Approaching 
the soliton core centered at x = Xg, the phase <p{x Xg it 0,t) matches the angles ±6, which 
delimit the chordus part of the trajectory, see FigUl From large scales we view only a jump 
if{x Xs,t) ^ 9{t)sgD.{x — Xg), which effect can be easily extracted if we generalize the scheme 



suggested earlier for static solitons at presence of interchain coupling 
phase mode is 
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25( 1 ■ The action for the 



Ssnd[f{x, t), e{t)] = ^ 1 1 dxdt {{dt^/uf + {d^^ - 29{t)6{x - Xs)f) 



(7) 



where u is the CDW phase velocity: u = gphCVph/v 4:Tt. Integrating out (p{x,t) from 
exp{—Ssnd['&L§]} ^ given 9{t), we arrive at the typical action for the problem of quantum 



dissipation 



Sd.s[e] ^"-011 dtdt' (^4^)' = -^-^11 'itdt'minit-tW) (8) 

that is Sdis ~ Yl I'^ll^o;!^- (The above equivalence of two forms for Sdis[G] holds for closed trajecto- 
ries of 0{t).) The dissipation comes from emission of phase phonons, while forming the long range 
tail in the course of the chordus soliton development. In real time, the dissipation comes from 
emitting the phase fronts (p{t,x) ~ Sgn{x)9{t — \x\/u). 



III. ONE- AND TWO- PARTICLE TUNNELING PROCESSES 

Consider the system of two chains (j = a, b) weakly coupled via hybridization of their electronic 
states, which amplitude is t± <C Aq. The chains are maintained at the electric potential difference 
U > 0. Electrons jump between the chain a at the potential U/2 and the chain b at the potential 
—U/2. Initially both chains are in their ground states with the same total number of particles 2M 
at one chain. In principle, this is a nonequilibrium situation with the total energy difference of 
2MU. By definition, we do not allow for an extensive 5M ~ M interchain charge transfer to form 
a common Fermi level: that would be an inhomogeneously reconstructed junction state [s], which 
requires for a special treatment. Considering rare events of one or two particle tunneling, we need 
to follow only the potential difference for their intragap states. 

For the interchain tunneling Hamiltonian 

H±= t± / dxi^l{x,t)^b{x,t) + ^l{x,t)^a{x,t)) (9) 

<a,b> '' 

the transverse electric current operator is 

i± = - ^l{x,t)-^a{x,t)). (10) 
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The average current (number of electrons per chain per unit length per unit time) is given by the 
trace of this operator over fermions and the functional integral over lattice deformations. 



j{U)= j j DAaix,t)DAbix,t)Trj^exp 



-SiAa) - SiAb 



dtH± 



(11) 



We shall need j{U) only in the lowest nonvanishing order of expansion over tj_. Each application 
of operators jj_,H± at a time t produces following effects: changing the occupation numbers 
of split-off states I'j by ±1, and multiplication by the product of corresponding wave functions 
^^{x,t)^{y{x,t), adding the voltage energy gain U. 



A. One electron tunneling. 

Consider the single particle channel in interchain tunneling for a standard ICDW, with a double 
g = 2 spin degeneracy of electronic states. 

Processes originated by the transfer of one electron between the chains appear already in the 
lowest second order of expansion of the current in powers of t±. The contribution Ji to the 
interchain current can be written as a functional integral over fields Aj, j = a, b: 

Ji-tlj d{x-y) j d{h-t2) J D[Aj][%{x,h)^hix,h)^a{y,t2)niy^t2) 

exp(-S'i(ti,t2, Aj(x,i))) 

were ipj are wave functions of the particle added and extracted in moments ti,t2 at fluctuational 
intragap levels Ej. Here the time dependent action S{ti,t2) describes the process of transferring 
one particle from the doubly occupied level -E^ < of the chains a to the unoccupied level E), > 
of the chain b at the time ti , and the inverse process at the time t2 ■ We have 



Si{h,t2,Aj{x,t)) = Sa{0,-OO,ti) + Sb{0,-OO,ti) + Sa{0,t2,oo) + Sb{0,t2,Oo) (12) 

+ Sa{-1, tl, ts) + Shih hM) + U{h- t2) (13) 

where Sj{i',ti,t2) = Sj{i',ti,t2; Aj) are the on-chain actions The branches' components are 
shown in the table below ( 6 stays for 0" or for tt — 0^), signs it correspond to a and 6): 



E/Ao 


Vo/Aoit<ti,t>t2) 


Vi/Ao ih<t< t2) 


-cose±U/2 


-^6lcos6l + ^sin6' 


(1 - fe*) cos + 1 sin 6* - C//2 
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FIG. 4: Branches for the one electron tunnehng. The potential evolves following the ascending solid line 
of Vb, then switching to the descending dashed line of Vi until it touches the axe ^ = Energy levels are 
shown as dotted lines. 

The corresponding branches are shown at the FigUl At this plot we have chosen a critical, 
minimal value of C/ = 2Ws = 1.27Ao (the thick line touched the zero level), that is 9i = arccos 2/7r ~ 
50°. 

The main contribution to the subgap current comes from extremal instanton trajectories, which 
minimize the action Si (1130 over both Aj{x,t) and r = t2 — ^i, at a given U, thus yielding 
minS* = S{U). The exact extremal trajectory is defined by equations 

6S/6Aj{x,t) = ; dS/dh = , dS/dt2 = 0. (14) 

The first equation in (fn|) becomes, in our projection, 6S/66j{x, t) = 0, hence the same equation 
for 6a and n — 6b- 



Itt 



The equation (jlSp describes the dissipative motion of an effective 6— particle with the variable 
mass 2/ in the effective potential — Vj^(t)(^), which depends explicitly on t via jumps of = u{t) at 
the tunneling moments ti. 

The second pair of equations in (I14p can be written, in our model, as 



Eoa{ti) + EQb{ti) = U ; EQ{ti) = Ao cos d{ti) = Aqcos^^ (16) 

This equation shows that the tunneling takes place when quantum fluctuations at both chains 
simultaneously create such potential wells, that they produce split-off intragap states satisfying 
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the resonance condition Ef, — Ea = U. (The last statement is Hterally vahd only within our model 
neglecting Coulomb interactions, see examples of these corrections in [l2. [isl.) 

Based on earlier results {ill , [isl augmented by reasonable assumptions (see Appendix), we 
conclude that, up to pre-exponential factors, the tunneling current is pro portional to the square of 
the single particle spectral density, which has been already calculated jlll. Iisl. as the PES intensity 
IpES- Ji oc t\llj^g{VL = U/2). E.g. near the threshold U 2Was 

U-2WasV^ ^ VF 



.m) = Miy^^^ .^=^»i (17) 

This law gives the vanishing current at the threshold, which is different from low symmetry cases 
[l^ . [27I where the cutoff was finite. A physical origin of this suppression comes from efficient 



emittance of gapless phonons in the course of tunneling, which drives dynamics to the regime of 
quantum dissipation [ill, [3] • 

Calculating the prefactor A requires for integration over Aj(x,t) around the extremal, taking 
into account the zero modes related with translations of the instanton centers positions Xj{t) and 
with the phases (pj = ip{Xj). The prefactor A = A{U) in ()17p is a slow function of U; it may 
yield powers of {U — 2Was)~^ with ^ ^ 1 [3 which dependence is negligible in comparison with 
the big power 2(3. But the U dependence of the prefactor becomes important in the 3D regime 
near the threshold. The 3D long range ordering, which energy scale is measured by the transition 
temperature <^ Aq, returns the tunneling to the normal dynamics. Then the peak of Ji{U), 
instead of zero (jlTp . will develop in a narrow vicinity < U — 2Ws < Tc of the threshold, in a 
qualitative agreement with experiment Q]- The factor {Tc/Wg)'"^'^^ will give the overall reduction 
of the current in comparison with the one expected for free electrons. 

Behavior near the free edge ~ Aq is governed by small fluctuations i] of the gap amplitude 
|A| = Aq + r] and of the Fermi level SEp = '~p'vpl2 via the phase gradient i-p' = dx'-p- The gap 
fluctuations dominate over the phase ones, so the results are similar to the ones for commensurate 
CDWs. A thorough discussion was given in |l3|. 



B. Spinless case. 



The spinless case, the degeneracy fif = 1, applies to spin polarized CDWs in high magnetic fields, 
which is a current experimental trend [26] . It describes also the XY model of the spin-Peierls state 
in magnetic chains. Also this case provides a tutorial to the regime of full phase slips at small 
[/, which otherwise appear only in higher forth order in t_|_ (see the next section). The interchain 
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jump of the fermion takes place between the filled level rising from E = — Aq and the empty level 
descending from E = Aq. The branches' components are shown in the following table and drawn 
at the Figl5] (all energies are in units of Aq). 



E 


Vo {t<h,t>t2) 


Vl {h<t< t2) 


-cos0±U/2 


-i6'cos6i+ isine* 


(l-i6i)cos6' + isin6'-J7/2 




FIG. 5: Branches for the spinless case. Here U = 0.2, then 6i — arccos(O.l) = lATrad. Initial branch Vq 
a,t 9 < Oi is given by the solid line, the branch Vi &i 9 > 9i - by the dashed line, energies are shown as 
doted lines. 

At [/ ^ the tunneling point moves to 9i vr/2, where the configuration is the amplitude 
soliton. It provides the exactly midgap state, which is necessary at U — > 0, hence there is no low 
boundary for U. We arrive at the result 

(The exponent is twice smaller in comparison with (jl7p because coefficients of each term in the 
action for g = 1 are 1/2 of those for g = 2.) 
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C. Bi-electronic tunneling. 



The joint self-trapping of two electrons allows to further gain the energy, resulting in stable 
states different from independent ASs: there are the phase slips leading to phase 27r solitons. An 
advantage of tunneling experiments is the possibility to see the phase slips directly, at voltages U 
below the two AS threshold U < 2Was of the pair-breaking, i.e. within the true single particle 
gap. The probability of the bi-electronic tunneling is relatively small as it appears only in the 
higher order ~ in interchain coupling. But it can be seen as extending below the one-electron 
threshold, where no other excitations can contribute to the tunneling current. 

Consider the forth order (~ t^) contribution to the current due to the possibility of electron 
pairs tunneling between chains. 

4 

J2^t\j n DAj{x,t) j Y[dyidtiexY>{-S2) 

j=a,b i=l 

[%{yi,ti)^a{yi,tiM{y2,t2)^a{y2, t2)K{y3, h)^b{y3, t3)K{yA, U)^b{yA, U) 

-{y4,t4< — > 2/3,^3}] (19) 



5*2 = 5'2(tl..4, Aj (x, t)) = Sa{0, -00, ti) + ^^(O, -00, ti) + Sa{0, t^, Oo) + S'fe(0, t4, Oo) 
+Sa{-1, tut2) + Sb{l, tiM) + Sa{-1, hM) + Sb{l, h, U) - U{t2 - tl + U - t^) 

+Sa{-2, tiM) + Sb{2, ti,t2) - 2U{t2 - h) (20) 

The action (pOj) describes the following sequence: The first jump takes place at t = ti, 9{ti) = 61 
from the twice filled level rising from E = —Aq to the empty level descending from E = Aq. Its 
branches are the same as shown above in the table for the one-electron tunneling. 

The second jump takes place at t = t2, 0{t2) = 62 from the singly filled level of the chain a to 
the singly filled level of the chain b. The branches' components are shown in the following table 
and plotted at the FigEl (All energies are in units of Aq.) 



Vo (t <ti,t> U) 


Vl {ti<t< t2),{U <t<ti) 


V2 {t2<t< h) 


-^9cos9+ ^sinO 

7T TT 


(1 - |6i) cos6'+ f sin6i- C//2 


(2- |6l) cos6'+ |sin6l- [/ 



Rigorously, the extremal trajectory is defined as above in Sec illl A] from equations 

5S , dS 

__ — _ = 0,j=a,6; — =0,^ = 1..4 

O0j{x,t) Oti 
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FIG. 6: Branches for the bi-electron tunneling. Here U = 0.4, then the transition (lines crossing) takes 
place at 9 — arccos(0.2) = 1.37 - below the amplitude soliton value 9 ~ 7r/2. Split-off energies are shown by 
dotted lines. The potential evolves following the thick lines: solid, then dashed ones. The transition of each 
of two electrons proceed via the intermediate branch Vi - the thin line. 

Similar considerations argue that equations are consistent provided 

Oaix, t) = 7T- Obix, t) ; Eo{ti) = tl = t2 , h = U (21) 

The last condition ti = t2 , = tells that only processes of simultaneous tunneling of pairs 
of particles contribute to the extremal action (j20p . The branch with = 1 is strongly virtual, it 
is passed over the microscopic time Co/u, hence at the FigElthe interval of Vi is reduced to the 
point. The intermediate state of = 1 gives no contribution because of the accidental degeneracy 
of our model. This is the artifact of the Peierls model: due to the conditions (j2ip the curves Vq, 
Vi — U /2, V — U intersect at the one point 9 where Eq = U /2. At C/ ^ the tunneling points 
move to 9i^2 '?i"/2, where the configuration is the amplitude soliton. 

Finally we arrive at the power law with the index vp/u twice that for the single particle process 

J2 oc exp [- 25] oc f — j (22) 

The bi-electronic channel in ICDW stretches the PG down to the small energy tx Tc, i.e. even 
to U = in the ID limit. The reason is that here the particles with the charge 2e - the phase 
2tt solitons are neither bipolarons, nor pairs of kinks as for systems with lower symmetries, see 
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[s, |9|. The tunneling takes a form of a coincidence of opposite ±2tt phase shps taking place 
simultaneously at adjacent chains. In the regime with the 2D or 3D long range order, the current 
J2 will vanish below the low threshold 2Wps where the characteristic intensity will be of the order 
of Jaocti (r,/Aor/". 



E 

O 

> 

■D 



2A peak 



-200 



-100 



Amplitude Solitons 




-50K 
55K 



100 



200 



V(mV) 



FIG. 7: (Color on line.) Experimental data from the collection of Latyshev at al, e.g. [gl, and their 
interpretation. The features correspond to the second, low T, CDW in NbSea. They are taken just below 
the transition temperature {T — - green curve), and at slightly lower T = 50K - red curve. The 

evolution of the AS feature from an edge to a peak corresponds to our guesses on the role of the long range 
order. 



IV. CONCLUSIONS 

We conclude that the three regimes are expected theoretically for the internal interchain tun- 
neling in ICDWs. 

1. High energy single-particle peak at 2Aq corresponding to the electron- hole excitation across the 
rigid CDW gap. 

2. Middle energy edge at 2Was ~ 4/3Ao corresponding to creation of two Amplitude Solitons, 
which are, unlike undressed electrons of the case 7^1, the exact pair-breaking eigenstates of the 
whole system. 

3. Low energy edge at Wps ~ Tc, which drops down even to U = in the ID regime, when the 
scale of the interchain ordering energy Tc may be neglected. 

The regimes ^^1 and #2 can be accessible also in optical absorption, while the regime ^^3 is 
specific to tunneling. Some more studies are necessary to reach a quantitative comparison with 
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experiment. There are effects of a finite temperature and of direct electronic interactions, calcu- 
lating the prefactors, taking into account the inhomogeneous distribution of the electric potential 
over a junction. Nevertheless, the general classification, assignment of features, and their location 
are already well covered by the theory, cf. P, Q]. These relations of theory and experiment are 
illustrated at the FiglT] 
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V. APPENDIX 

A. Selftrapping details. 

In the following the units of energy and length are Aq ^ 1 and = ft-vp/^o ^ 1- 
Topologically nontrivial trajectories in the plane of the complex order parameter A are conve- 
niently parameterized by a family of chordus solitons 



Asoi = Eq- iko tanh(A;ox) (23) 

Here the split-off energy Eq and the inverse localization length ko are parameterized by the chiral 
angle 29 (see Fig(T| as 

Eq = ziz COS 9 , ko = sin 9 ; < 9 < tt 

with signs it corresponding to two directions of the trajectory (the positive one is shown at Fig{T]). 
We shall assume that vr > ^ > 0, so that the negative chirality is obtained by transformation 
9 ^ TT — 9. The level Eq is filled hy < u < g fermions, while vacuum levels at E < — Aq are all 
filled with the total spin degeneracy factor g. (In applications, g = 2 - electrons, ox g = \ - spinless 
fermions.) The spectral flow provided by the positive phase slip draws the energy level Eq = cos 9 
down from Eq = 1 at 9 = towards Eq = —1 at 9 = it. Its energy branch is 

V^{9) = {u - -9)cos9 + -sm9 , Eq = cos9 , u = 0, ...,g (24) 
TT vr 

The spectral flow provided by a negative phase slip draws the energy level Eq = — cos 9 up from 
Eq = —1 at 9 = towards Eq = 1 at 9 = tt. Its energy branch is V-u{9) = V^{tt — 9). 
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FIG. 8: Schematic plots of time dependencies for 9 and 6 for one-electron tunneling. 



B. Dissipative dynamics. 



Here we shall analyze general properties of the extremal solutions. The total on-chain action 
can be written as 

vp/u 



s{e,y) = j dt[f{e)e^ + v,{e)] 

Its variation, Eq. ljlSp . can be written as 



27r2 



dtdt'e{t)\n\t-t'\9{t') 



(25) 



2f9 + -V + 



-oo (t - t') 







(26) 



Here 9 means the time derivative, while ' means the derivative over 9. 

Consider the single-particle process when jumps from to ±1 and back at times ti^2- The 
expected time dependencies for 9 and 9 are shown schematically at the FiglHl 

Let r = t2 — ti — > oo, then the solution should decouple in two solitons evolving between 
=^ ^0 ^-iid ^0 =^ 0) s-iid we may consider only the the first one. Here 9q is the stable exited state 
corresponding to the minimum of Vi : 9q = ir/l for the spin 1/2 case and 9q = n for the spinless 
case. Contrary to usual problems with degenerate vacuums, here there is a discontinuity at the 
floating point ti when < < ^o- At the level 9i, V has a jump, hence the jump in 9 and the 
cusp in 9. The result will confirm that the total increment of 9: A9 = 9q is accumulated in a finite 
vicinity of ii, then Eq. (j26p becomes 

VF 9o 



2f9 + f'9^ - V 



(27) 



vr^n t — ti 

Here 9 means the time derivative, while ' means the derivative over 9. At the left asymptotics 
t — ti — > — oo, V = Vq ^ 9'^g/2>'K and f ^ 9, then the l.h.s. of (j27p is dominated by the term V' and 

— 1/2 

we obtain the asymptotics 9 ^ \t — ti\ ' . The right asymptotics at t — ti — > oo depends on the 
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degeneracy. For g = 1, the final state 9o = it is the full phase slip, equivalent to the initial state 
6 = 0; hence the right asymptotics is the same as the left one. For g = 2, 9q = Tr/2 corresponds 
to the amplitude soliton. Here / cnst and Vi ~ ^ + ^ {0 — ^^r)^; again the l.h.s. of (j27p is 
dominated by the term V and we obtain the asymptotics 6 — Oq \t — and similar near t2- 

The obtained asymptotics show that the integral in the dissipative term of Eq. (j26|) is convergent 
around the impact points tj, which makes our assumptions to be consistent. 

For a finite r we have a two-soliton process, and, in the leading order the solution is the sum 
of two contributions. Then for the outer interval |t| ^ r, 6^ ~ rt-3/2. Within the inner interval 
ti <t <t2 

We find for the turning point Omx — do ~ ^/t- 

The same procedure can be applied to calculating the action ([25]): again, the integral is concen- 
trated at t, t' being near the points ti, t2- When both t, t' belong to the same vicinity, we obtain a 
constant contribution to the action which may be omitted. When t and t' belong to vicinities of 
different points ti or t2, then the contribution is ~ Inr. Collecting only the contributions which 
increase with r, we obtain the action St, its extremum Sextr = S{Textr), and finally the tunneling 
probability Ji ~ exp(—2Sextr)- 



1 I ' 1 ' extr 1 I ,j-r J J \ ) 

U \TT J U \-K J {U - Umin) 



U — Umin \ n Vp f Oq^"^ 



V Umin J M V vr y 

Consider finally the question of uniqueness of the chosen solution. Phases 9j obey (up to the 
opposite rotation 6 ^ tt — 9) same equations, characterized by common switching moments ti. 
Their solutions could be different in principle, as characterized by different sets of switching phases 

— ^"'(^i)) * = 1,2 and j = a,b. Eqs. (jl6p reduce twice the number of boundary variables, so 
that e.g. only 9f are independent unknowns. The instanton process starts at t = — oo with values 
9^ = 0, comes to the moment ti at values 9f and 9\ related by ()16p . reaches the turning points 9mx 
at turning moments tmx', then it comes to the moment t2 with values ^2 and ^2 related by (|16p . 
and finally returns to values ^-^ = at t = +00, see FigUl We shall assume that the process is 
symmetric, as it is demonstrated by solutions of nondissipative models, see {l^. [isl. and of simpler 



purely dissipative ones 
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33]. Then 6*^ = 6ij, and t^x = (^i + t2)/2 becomes the same for 
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both chains. We are left with a single unknown boundary phase at each chain, uniquely related to 
time interval t = t2 — ti. The parameter is common for both chains, then the solutions 9^{t) are 
identical. In view of (|16p . they will be determined by the same condition £'oa(^j) = EQf,{ti) = U/2. 
We arrive at the effective one chain problem with a doubled effective action for the current. 
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